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(Dated: February 7, 2008) 

Abstract 

We report on the observation of the rf c (2 1 So), the radial excitation of ?? c (l 1 So) ground state of 
charmonium, in the two-photon fusion reaction 77 — > rf c — > K^K^ir^ in 13.4 fb _1 of CLEO II/II.V 
data and 9.2 fb _1 of CLEO III data. The data have been analyzed to extract the rf c resonance 
parameters. 



* Submitted to the International Europhysics Conference on High Energy Physics, July 2003, Aachen 
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I. INTRODUCTION 



Knowledge of the hyperfine (spin singlet-triplet) splitting is important for the under- 
standing of the spin-spin interaction in quarkonia. The existing experimental knowledge is 
limited to the 117±2 MeV splitting between J/ip(l 3 Si) and ^(l 1 ^) It is important to 
find out the magnitude of the corresponding splitting between the radial excitations tp'(2 3 Si) 
and rj' c (2 1 So), which samples more confinement-dominated part of the cc potential whose spin 
dependence is open to speculation. The proximity of if)' and r]' c to the DD breakup threshold 
also brings in channel-coupling effects. Several searches for rj' c have been reported in the 
literatureQ-jij]. 

The Crystal Ball Collaboration at SLAC reported the identification of 7]' c in the radiative 
decay of ^'H- Structure due to a purported 7 of energy 91±5 MeV was interpreted as 
the signal for M(r]' c ) = 3594± 5 MeV, and a 95% confidence level upper limit on its total 
width was set as 8 MeV. The E760 experiment at Fermilab made an unsuccessful search for 
r)' c in the reaction pr> —> rj' c —> 77 in the mass region 3591-3621 MeV with ~6 pb _1 of 
invested luminosity [3|. They set a 90% confidence upper limit for B(rf c — > pp) X T 77 (^) of 
~0.7 eV for assumed total widths between 5 and 10 MeV. The search was repeated by the 
successor experiment E835 with the measurement with ~30 pb _1 of luminosity. No evidence 
for r]' c was found in the mass region 3575-3600 MeV, and a 90% confidence upper limit for 
B(r)' c — > pp) x T 77 (^) of ~0.7 eV was again set for assumed total widths <15 MeVjiJ]. It is 
worth noting that the corresponding value for rj c is B{r\ c — > pp) x r 77 (?7 c ) ^9 eV, i.e., about 
an order of magnitude larger. 

The next search for rj' c was done by DELPHI in the 77 fusion reaction e + e~ — > 
e+e"(77) e + e~(r]' c ), rf e -> p°p°, K° s Kn, K*°Kn, K° s K s nn and K + K~K + K~ . No sig- 
nal was observed in the mass range of 3500-3800 MeV, and a 90% confidence upper limit 
r 77 (?^.)/r 77 (77 c ) < 0.34 was set||| assuming equal decay branching ratios for rj c and T]' c . In a 
similar search L3 established a 95% confidence limit of r 77 (^)/r 77 (?7 c ) < 0.29j6j. 

Recently, the BELLE Collaboration has reported an rj' c signal with >6a significance in 
the decay of B — * Kr)' c , with v( c — > KgK ± ir T , from which they determined M(f]' c ) = 
3654±6(stat)±8( y syst) MeV0, and T tot (r]' c )= 15 t^(stat) MeV. The BELLE Collaboration 
also presented evidence at the 3.4(7 level for an r)' c signal in the reaction e + e~ — > (J/if))r}' c , 
with M(i c ) = 3622 ± \2{stat) MeV@. 

Theoretical predictions for the mass and width of the T]' c meson have been so far based on 
potential model calculations and an analogy between the radial excitation pairs (J/i/j)/i/j' 
and T] c /r]' c . The predictions for masses range from M(j]' c ) = 3594 MeV to M(j]' c ) = 3629 
MeV, and the predictions for the ratio of the two-photon partial widths, r 77 (^)/r 77 (r7 c ) 
range from 0.38 to 0.77. Thus, the current best experimental results for both the mass and 
two-photon partial width of the r)' c are at variance with theoretical predictions. This makes 
for great interest in determining these parameters with precision. 

At CLEO we have searched for the r)' c meson in the final state KgK ± 7c T via 77 fusion in 
e + e~ data at the T(4S) resonance and its vicinity, using 13.4 fb _1 of CLEO II/II.V data. 
Having observed strong evidence for the r)' c , we searched for its corroboration in 9.2 fb _1 of 
CLEO III data. The observation of r] c in the same reaction and decay channel was reported 
by CLEO earlier @. 

Preliminary results of the present study were presented earlier The BaBar Collabora- 
tion has also recently pre sented their preliminary results of the study of r\ c and rf c production 
in the same reaction 
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The layout of the paper is the following. Section II describes our data samples and event 
selection criteria. Section III discusses our analysis procedures, while Section IV presents 
our estimates of systematic uncertainties. A summary of the results is given in Section V. 



II. DATA SAMPLE AND EVENT SELECTIONS 

The CLEO II/II.V data sample consists of 13.4 fb -1 of e + e~ luminosity. Approximately 
one-third of the data were taken with the CLEO II configuration of the detector, and two- 
thirds with the CLEO II. V configuration. The total luminosity of CLEO III data used in 
the current analysis is 9.2 fb -1 , ~ 70% of the CLEO II/II.V luminosity. However, this newer 
data have the advantage of the superior hadron identification, as described below. 

The detector components most useful for this study were the concentric tracking devices 
for measurements of charged particles, operating in a 1.5 T superconducting solenoid. For 
CLEO II[12|, this tracking system consisted of a 6-layer straw tube chamber, a 10-layer 
precision drift chamber, and a 51-layer main drift chamber. The main drift chamber also 
provided measurements of the specific ionization loss, dE/dx, used for particle identification. 
For CLEO II.V . th e innermost chamber was removed and replaced with a three-layer silicon 
vertex detector |l3j]; in addition the gas in the main drift chamber was changed from a 50-50 
mixture of argon-ethane to 60-40 helium-propane. 

For CLEO III Hi, the entire charged particle tracking system was replaced. The new 



tracking system consists of four layers of double-sided silicon strip detector, surrounded by a 
new, 47-layer drift chamber with a smaller outer diameter (l5|. Careful design optimization 
led to the momentum resolution of this tracking system comparable to that of CLEOII/II.V. 
The CLEOII time-of-flight system was eliminated and in the s pac e created by it and the 
smaller drift chamber a ring imaging Cerenkov detector (RICH) [16] was installed. 

The RICH distinguishes charged kaons from pions over 80% of the solid angle, which is 
smaller than hadron identification with dE/dx. For charged tracks with momenta below 2 
GeV/c (the momentum range of most of our kaon candidates) the RICH identifies kaons with 
efficiency greater than 81% while having less than 2% probability that a pion is misidentified 
as a kaon. 

In addition, the trigger for CLEOIII was improved in its flexibility, efficiency, and redun- 
dancy. Given that it relies on fewer types of devices with better design than in CLEOII, the 
trigger is also easier to reliably simulate. 

The Monte Carlo simulation of CLEO detector response was based upon GEANT. Simulated 
events were processed in the same manner as the data to determine the K s K ± ir Zf detection 
efficiency Efficiencies obtained in this fashion are listed for major selection criteria for 
both CLEO II/II.V and CLEO III in Table [I] The overall detection efficiency for rj' c is 
13.8%(CLEO II) and 11.6%(CLEO III). 

The data sample used contained four charged particles, with at least one K s candidate. 
The Kg vertex was reconstructed from its decay to 7t + 7t~ and was required to be displaced 
from the e + e~ interaction point. The mass of Kg candidate was also required to be within 
0.490 < M(nn) < 0.505 GeV/c 2 (0.488 to 0.508 GeV/c 2 for CLEO III). Furthermore, the 
Kg momentum vector was required to point back to the interaction point. The efficiency of 
these criteria, as determined from our simulations, was 77% for r] c and 79% for rf c . As shown 
in Table I, the CLEOIII efficiencies were somewhat smaller than those of CLEO II/II.V. 

The two charged particles other than those from K s decay were examined for the hy- 
potheses of being either kaons or pions. In CLEO II/II.V, this was done by using the dE/dx 
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TABLE I: Efficiencies for detecting the r\ c and rj' c signals in %. Efficiencies for asterisked cuts were 
obtained for events which have passed all other cuts. 





CLEO II/II.V 


CLEO III 




Vc 


Vc 


Vc 


Vc 


4 reconstructed tracks 


35 


40 


34 


37 


Ks identification 


77 


79 


61 


64 


K/tt particle ID * 


90 


93 


75 


78 


P T < 0.6 GeV/c * 


93 


95 


84 


84 


Small E neut * 


84 


83 


82 


80 


Trigger * 


69 


73 


99 


99 


Overall 


10.0 


13.8 


9.8 


11.6 



and TOF information (when available), and by applying appropriate probability criteria. 
The efficiency of this particle identification, again determined from our simulations, was 
90% for the rj c and 93% for the r)' c . 

For CLEO III, all events were used in which the charged kaon candidate is identified as a 
K by the RICH detector. When the K candidate was not identifiable by the RICH (mostly 
when it is outside the RICH fiducial volume), it was identified by dE/dx measurements from 
the drift chamber. When its momentum lay between 1 and 2 GeV/c, however, it was difficult 
to distinguish K from a charged pion using dE/dx and we did not use those K candidates. 

Other criteria used included using specific trigger requirements, requiring a small trans- 
verse momentum of the KgK ± n T system (P^<0.6 GeV/c), and requiring the sum of neu- 
tral energy not associated with the charged particles be minimal, (E neut <0.2 GeV for 
CLEOII/II.V; <0.3 GeV for CLEO III). 

III. ANALYSIS 

The KgK^ir^ invariant mass plot using our final event selection for CLEO II/II.V is 
shown in Fig. la. Enhancements at masses ~2985 and ~3643 MeV are visible. Henceforth 
we call these r] c and T]' c enhancements, respectively. 

In order to extract numerical results from these data we have made maximum likelihood 
fit to this spectrum using a second order polynomial background and two Breit-Wigners, 
with their parameters kept free. This fit is overlaid in the Figure. As presented in Table HTl 
the fit gives 36+11 events at a mass of M(r)' c ) = (3642.7±4.1(stat.)) MeV. A fit using a Breit- 
Wiegner convolved with a double Gaussian representing the resolution function extracted 
from Monte Carlo for each of the signals is in good agreement with the fit using Breit- 
Wigners only. Note that we assumed that the potential interference between the rjciv'c) 
signal and the continuum two-photon production of KsKn is negligible when we equate the 
peaks of the distribution with the masses of the r)V\ The same assumption was made for 
the rate measurements. 

The "significance levels" of the enhancements in this Table were obtained as 
a =\J — 2ln(L Q / L max ), where L max is the maximum likelihood returned by the fit described 
above, and Lq is the likelihood returned by the fit with either no rj c or no rj' c resonance. The 
search for rj' c resonance was made in the mass region 3580-3700 MeV. Since the mass and 
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M(K s Kjr)(GeV) 



FIG. 1: KgK^Tr^ invariant mass in the reaction 77 — > KgK^ir^ from the (a) the CLEO II/II.V 
data and (b) the CLEO III data. The curves in the figures are fit results using a second order 
polynomial for the background. For CLEO II/II.V two Breit-Wigners are used for the r\ c and r]' c 
enhancements; for CLEO III the fit shown uses Breit-Wigners convolved with Gaussians. 

width are not a priori known, a more appropriate way of estimating the signal significance 
is to determine the probability that the counts for the background alone could statistically 
fluctuate to the level of the counts of signal plus background. It is found that the significance 
levels determined in this manner is smaller than those listed in Table ITT1 bv about 0.5a. 

Fig. 1(b) shows the preliminary K^K^Tt^ invariant mass distribution from CLEO III 
after all criteria were applied. It shows clear enhancements at the r\ c mass as well as at the 
mass at which the CLEO II/II.V data indicated the r( c . This is strong confirmation that the 
rj' c signal we observed in the CLEO II/II.V data is real and not due to statistical fluctuation 
or some other artifact. 

We fit the mass spectrum with second order polynomial representing the background 
under the rf c candidate peak and a Breit-Wigner function convolved with a double Gaussian 
representing the instrumental resolution of 5.7 MeV (96%) and 40 MeV (4%) as determined 
from our Monte Carlo simulation. We have also used power-law background and exponential 
background shapes and averaging over all three background shapes, we obtain preliminary 
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TABLE II: The results of the analysis. The yields of events and the masses are from the fits to the 
data. The significance of the enhancements are based on changes of likelihood as described in the 
text. The CLEO III results are preliminary. 





CLEO II/II.V 


CLEO III (preliminary) 




Vc 




Vc 


rfc 


Yield (events) 
Mass (MeV) 
significance 


287 ± 28 
2984.7 ± 2.1 
15.9cr 


36 

3642.7 ± 4.1 
5.0a 


203±22 
2982±2 
14.2cr 


29±8 
3642.5±3.6 
5.7a 



results of (29±8) events in the peak at a mass of (3642.5±3.6(stat)) MeV. These fit results 
are also shown in Table |HJ This measurement of the vj c mass is completely consistent with 
that obtained from the CLEO II/II.V data. The statistical significance of the rj' c signal is 
determined to be 5.7a, using the likelihood method described above. Note that in this case, 
this significance is the appropriate quantity since we search for the r)' c at a particular mass, 
namely that determined in the CLEO II/II.V analysis. 

Photon-photon fusion is expected to populate positive charge conjugation resonances 
mainly when the photons are almost real, i.e., when the transverse momenta of both of 
them, and therefore of the sum of final state particles is small. In order to examine whether 
the observed rj' c peaks are due primarily to two-photon events, we examined the production 
of r)' c in two subregions of transverse momentum: < P T < 0.2 GeV/c and 0.2 < P T < 0.6 
GeV/c. Table ITTT1 shows the results of this study. The CLEO II/II.V and CLEO III results are 
statistically consistent with the expectations from our two-photon Monte-Carlo simulations. 

TABLE III: The transverse momentum dependence of production. 



Pt 


0-0.2 GeV/c 


0.2-0.6 GeV/c 


2 7 Expectations 


74% 


26% 


CLEO II/II.V 


(55 ± 18)% 


(45 ± 18)% 


CLEO III 


(62 ± 19)% 


(38 ± 19)% 



Also of interest is the two-photon partial width of the rf c as compared to that of the r\ c . 
The ratio of the number of events of r)' c and r\ c is: 

NijQ - hMl x B Wc -> gggjO x $ ("Vc) x <v'c) (1) 

N(r) c ) r 77 (^ c ) x B( Vc K s Kir) $(m„J e{r) c ) ' 1 ] 

Here, for the CLEO II/II.V analysis: 

N(r)' c )/N(r] c ) = 0.125lo;o4i is the ratio of r]' c and r\ c events from Table ILT1 
$(m 7? /)/<I>(mj ?c ) = 0.42 is the ratio of the two-photon fluxes at the rf c and r\ c masses, 

as calculated using the known two-photon flux function (lj. This ratio has a systematic 

uncertainty of ±2%. 

e(j]' c ) / e(j] c ) = 1.38 is the ratio of efficiencies as calculated from Monte Carlo simulation, 
and is shown in Table HI 



7 



Using these values in the equation above, we obtain 



r 77 (r] c ) x B{r] c -* K s Kn) 

Performing this same calculation for the CLEO III data, we obtain R(rj' c /rj c ) = 0.29 ± 
0.09(stat), which is preliminary but consistent with the CLEO II/II.V result. 



IV. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties in the CLEO III results are still being investigated; therefore 
we only discuss those for the CLEO II/II.V analysis below. 



A. Mass of the rj c meson 

We estimate systematic uncertainties on our mass determinations with the following 
considerations. 

• It is difficult to estimate the contributions to systematic biases which arise from different 
choices of event selection when statistical uncertainties dominate. However, the statistical 
errors in the r\ c peak are much smaller than those of the rj' c peak, and we consider the extreme 
variation in r\ c mass in the different event selections we have attempted as providing an upper 
limit of 1 MeV as the systematic uncertainty due to event selections. 

• We have made a Monte Carlo study of how the masses of r] c and i]' c shift between the 
input values and those obtained after the detector simulation. For an input mass and width 
of M=2980 MeV, T tot =27 MeV, the shift, AM(r) c ), is 0.9±0.6 MeV. With an input mass 
and width of M=3640 MeV, r io4 =15 MeV, the shift, AM(r^), is 0.7±0.3 MeV. We conclude 
that the systematic error due to this source is ±1.0 MeV. 

• We have also tested how the masses of the particles that can be reconstructed compare 
with their known masses. We find that the reconstructed mass M (Kg)=A97 .69±0 .06 MeV. 
The PDGQ value is M( J FT°)=497.67±0.03 MeV. This small deviation, AM(K° S )=0.02±0.07 
MeV, is consistent with the conclusions reached for the mass of the J/ip meson, namely 
AM(J/^)=0.4±0.4 MeV, which was measured separately. We therefore conclude that the 
absolute calibration of our mass scale has an uncertainty of <1.0 MeV. 

• It is possible that the choice of background shape effects the mass determination. From 
the fits made with polynomial, power-law, and exponential backgrounds we find that M(r) c ) 
varies by <0.7 MeV, and M(r] f c ) varies by <0.1 MeV. We therefore estimate that this source 
contributes < 1.0 MeV to the systematic uncertainty in masses of rj c and rj c . 

Adding the above contributions in quadrature we determine the systematic uncertainty in 
the mass of rf c to be <2 MeV. We note, however, that r\ c mass obtained in the present study 
(Table II) differs from our earlier published result by ~4 MeV. Provided this difference, we 
assign the systematic uncertainty in our measurement as ±4 MeV, at present. Studies are 
being made to understand the mass difference, and we expect that systematic uncertainties 
at the level of ±2 MeV will be achieved in the final results. 
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b. R{n'chc) 



This value of R cited earlier was obtained using second order polynomial function for the 
background. The major source of systematic uncertainty in the ratio R is found to be the 
choice of background shape. The values in case of power-law and exponential backgrounds 
is essentially the same, R = 0.13^^l(stat) and R = 0.14^^1 (stat). As our final result for 
R we quote the average of the extreme values and assign the difference as systematic error. 
The systematic error due to efficiency ratio was ~2%, which is negligible. 

Thus, 



V. SUMMARY 



We have analyzed 13.4 ft)" 1 of e+e™ data of CLEO II/II.V and 9.2 ftr 1 data of CLEO 
III at T(4S) and its vicinity for resonances in the reaction e + e~ — > e + e~77 — > e + e~r]' c — > 
e+e-iK^K^). 

In the K° s K ±r K T invariant mass plot resulting from the CLEO II/II.V data, we see an 
excess of events at mass of ~3643 MeV with the significance level of over 4a in addition to 
i] c . We attribute this excess to the excitation of the rf c resonance. Observation of an excess 
at the same mass in the CLEO III data confirms the discovery of r]' c in the CLEO II/II.V 
data. 

Our results for the r)' c meson are as follows. 
The value of its mass has been found to be: 

M(^)cleoii = (3642.7 ± 4.1{stat) ± 4.0{syst)) MeV , (4) 

^"(^c)cLEOIII = 

(3642.5 ± 3.6(stat)) MeV(preliminary) , (5) 

thus, 

AM = M(if>') - M(rf c ) = (43 ± 4(stat) ± 4(syst)) MeV . (6) 
We also determine the ratio 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 
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